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APPEEDICES 



SIITOPSIS 


In this thesis problem of capacity expansion in 
single ccmnorlity and nulticonnod ity networks is considered. 
The problem considorod is to increase the capacity of 
network to o.cconiiiodato increased flow by increasing the 
capacity of existing arcs or by introducing new arcs. 

It is assunod that only cost under consideration is the 
cost of increasing the capacity of an arc. 

For single commodity flow network four heuristics 
and branch and bound method are developed. Computational 
porformancG of thoso algorithms is tested hj randomly 
generated problems. Computa-tional results of these test 
problems on IBM 7044/1401 system are presented. 

A branch and bound method for multicommodity flow 
network capacity expansion problem is developed. 



CHAPTBR I 


IBTRODUCTIOl 

Physical distribution systems such as clGCtricity 
or water supply, telephone or microwave links road or 
rail netw'orks , etc. can bo ne-turally considered as networks 
of flow. A common feature of these systems is that 
expansion may be necessary with lapse of time as demand 
requirements may change. In the case of an electricity 
power network, let us consider power flow from a specific 
source node ITg to a sink node and assume that the power 
requirement of the sink is increased from some value P^ to 
another value P2» The existing network may not be able to 
accommodate this increased power flow from to and 
what may be required is to add extra arcs to the network 
to make such a flow possible. Arcs can be added to the 
network at a specific cost and capacity, and the problem 
is then one of choosing which arcs to add so as to increase 
the maximum flow capacity (from Fg to ) to a value at 
least as .large as P2 and at the same time incur the least 
cost. 

If in a network there is only one source node and 
one sink node and flow is also of one kind only then this 
problem is referred to as single commodity flow problem. 
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If there are many sources and many sinks in the 
network and if we make the restriction that the flow from 
certain sources must be sent to certain sinks, then the 
problem is a multi -commodity flow problem. In case of 
single commodity flow problem, an undirected arc can always 
be thought of as two directed arcs of opposite direction, 
and arc flows of opposite directions can be cancelled. 

But arc flows of different commodities can not cancel each 
other. This is one of the central difficulties of the 
multi-commodity flow problem. 

In Chapter II,- we present literature survey for 
capacity expansion problem and other related problems like 
minimal cost flow problem and fixed cloarge problem. 

Chapter III deals with some of the pure network problems, 
which arise as subproblems in either one of the heuristic 
methods or the exact method, which are developed in this 
thesis. Some algorithms available from literature which 
were used in this work for these problems are also pre- 
sented. Computer codes for these algorithms were developed 
and are presented in Appendix B. 

In Chapter IV we present mathematical programming 
formulation and solution methodologies for single commodity 
flow network capacity expansion problem. The solution 
methodologies which were developed for this problem 



include four heuristic algorithms and an exact method. 

Large number of problems were solved. Computational 
results are also presented. 

Multi-commodity flow network capacity expansion 
problem is discussed in Chapter V. Mathematical programming 
formulation for the problem is presented. An exact method 
for the solution of the problem is suggested. 



CHAPTER II 


IITEPATURE REVIEW 

In this chapter we present review of the available 
literature for following problems, 

1. Shortest path problem 

2. Maximal floxf problem 

3. Minimal cost floxir problem 

4. Hotwork capacity expansion problem 

5. Multi-commodity flow problem, and 

6. Fixed charge problem. 

These problems are chosen for review of literature 
because few of them resemble closely with our problem and 
others arise as subproblems in various methods developed 
in this work. 

2.1 SHORTEST PATH PROBLEMS 

Consider a network (V, E), where V is the set of 

nodes and E the sot of arcs. Every arc A. - ^ E has asso- 

ij 

ciated with it a distance c.^. The problem is to find a 
path from any node to any other node with sum of the 

distances c^ ■ of arcs in the' chain a minimum. 
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More than 400 papers have been published in the 
area of shortest path problems and allied areas [43] and 
we have got the least intention to discuss them here. ¥e 
refer for a detailed bibliography on these to Pierce [43]. 

Some more recent works published in this area are 
those of Hadlock [23] and Fredma,n [15]. 

Hadlock [23] has developed a shortest path algorithm 
for grid graphs. Grid graphs arc a simple class of planer 
graphs for which the vertices can be assigned integer 
coordinates so that neighbours agree in one coordinate and 
differ by one in the other coordinate. Grid graphs arise 
in applications fi’om the layout design of integrated circuits 
to idealized models of city street networks. In many 
applications, a shortest path between two vertices is 
needed. The best known algorithms for the shortest path 
in a general graph of n vertices are of complexity 0 (n ). 
Hadlock *s algorithm however taking advantage of the concept 
of direction present in grid graphs is 0(n) in the worst 
case and Oifn) in the best case. 

2 

Fredman [15] has shown that ) comparisons and 

additions suffice to solve the all pairs shortest path 
problem for directed graphs on ^ vertices with non-negative 
edge weights. This result is exploited to produce an OCn-^) 
algorithm for solving the shortest path problem. 
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2.2 FLO¥ PROBLEM; 


The maximal flow problem can simply be described 
as finding maximum possible flow from source node to 
sink node llj. without violating the capacity constraints 


X. ■:< b- . on arcs. 

-L J i J 

A labelling algorithm (Ford a 
which was used to solve sub-problems 
described in detail in Section 3.1. 


jid Pulkorson [12]) 
in this thesis is 


A number of resoarchors Imve published their work 


in this area. Our intention will be to review only recent 
litera.turo published in this area. 

Dinic [11] has given an algorithm which solves the 
problem exactly in the general case after not more than 
On p (machine) operations, where n is the number of nodes 
of the net, p is the number of arcs in it and c is a 
constant not depending on the network. For integer data 
this algorithm, like the algorithm, of Ford and Fulkerson 
gives an integer solution with a supplementary estimate 
of the number of operations C^^np plus an estimate of the 
last. This algorithm is an iterative process like Ford 
and Fulkerson's [12] algorithm. 

ICarzanov [33] has presented an algorithm with an 
improved bound O(n^). Malhotra, Kumar and Maheshwari [38] 



7 


has given another algorithm, which is sitapler than 
Ivarzanov's [12] algorithm to solve the per stage flow 
problem in O(ii^) steps. 

2.3 MliTIML COST FLOW PROBLEMS: 


Consider a connected netTfork consisting of nodes F. 
and arcs leading from to . Among the nodes 
there is a special node called the source and a special 
node called the sink. The flow from to IT. in arc 


ij 


s . t . 


denoted 

by x.^.. 

ITow the minimal cost flow 

formulat ec3 2 ?:as i 


rain Z : 

= ^ c . . 

^ Id 

(x. . ) 

Id 



r -.-v for j = s 

^'x. . - 

^ Id 

k ^ 

J 0 for j ^ s, t 


pn 


(2.3.i; 


( 2 . 3.2 






V for j 


where aro non-nogative functions of x. and the 

arc flows are required to bo positive integers or zero. 
The parameter v, x/hicn should bo a, non negative integer, 
represents the total flow from source to sink. Fote that 
equations (2.3.2) express the conservation of flow at^^^ ^ 
nodes other than the source and the sink. The main 
computational procedures available for this problem whore 
linear functions of x^^- are the primal dual 



type algorithms of Ford and Fulkerson [14] and Busacker 
and Q-owen [4]. These are dual methods in which feasible 
flows become available when the computations terminate. 

Hu [29] has developed an algorithm where the cost 

functions c . - (x. .) are convene functions of x. .. In spirit 
1 J 1 J 13 

this aHgorithm is closely related to Beale's [2] and 
Busacker and G-owen's [4] algorithm. Hu's [29] algorithm 
is discussed in greater detail in Section 3.3 of this 
dissertation. In the later part of this paper Hu has shown 
that the problems associated with electrical networks with 
increasing capacity of a network under a fixed budget and 
with Max-flow Min-cut theorem may all bo formulated in to 
minimum cost flow problems in convex-cost networks. 

Fulkerson's 'out of kilter' algorithm (described 
in [14]) is essentially a primal method in that it can be 
started with a feasible flow or one becomes available at 
an early stage, 

Menon [40] has outlined a primal algorithm for the ; 
transhipment problem where transportation costs over each 
arc are convex function of amount shipped. 

Klein [35] has added one more primal method for both 
minimal cost flow and assignment transportation problems, 
with slight modification it can also be used for problems 
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invel-ving con-vex costs. It shares with other primal 
methods, the propertjr that it can he sta,rted with a ' good ’ 
solution and a hotter one is always a.vailahle in case 
early termination of computations is required. 

A more recent primal algorithm for minimal cost 
flow problems with convex costs is that of '¥Gtntra.uh [53]. 

2.4 NEB'JORIC CAPACITY BXPAYSIOIT PPOBIM ; 

LeBlanc [37] ii'i his paper has developed a long 
range planning model of transportation systems to assist 
planners in assessing the impact of va.rious levels of 
service in a transportation network. Es.i3. road distribution 
networks arc studied in particular, and emphasis is placed 
on the problem of determining optimsLl levels of service 
(improvements, degradations, abandonements ) on each arc 
in the existing network. The problem considered is to 
optimally design a transportation systemj the term design as 
used, does not coniprehond new networks, but more importantly 
the optimization or 'prunning' of existing transshijiment net 
works. The model is formulated so tha.t improvements to the 
shipping arcs are included as decision variables. The con- 
cave transhipment problem is extended to the case where 
coefficients of the, shipping cost functions are treated as 
decision variables. Thus the model determines optimum 
movements of frei^t between nodes and optimum improvements/ 
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degradations to the netwoi*k to minimize shipping costs 
plus na-intenanc e costs on the shipping arcs. The I'esult 
is a nonconveXj nonconcave transhipment prchlem which is 
then converted into a conco^ve minimization problem for 
which global- optimal solutions can be found. Unlike 
leBlanc [57], Christofides and Brooker [6] consider the 
problem of increasing the capacity of existing network 
where new arcs of finite capacity can be added at certain 
cost. The problon considered is s Given an existing network, 
a list of arcs which could be added to the network, the 
arc costs and capacities, and a.n available budget, choosing 
which arcs to add to the network in order to maximize the 
maximum flow from a source s to a. sink t subject to the 
budgetary constraint. 

The paper describes an efficient tree search algorithm 
using bounds co-lculatod by a dyneanic prograniaiing procedure 
which arc very effective in. limiting the solution space 
explicitcly searched. Computational results for a number 
of mediuiM sized probloms arc described and computing times 
are seen to bo reasonable. 

Howard and Homhauser [27] has considered this problem: 
given a sequence of predicted demands for 1 time periods, 
determine the optimal investment decision in each period 
to minimize a linear investment cost and a strictly convex 
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cost of capacitor . Tho relationship hotweon capacity and 
the investment decisions is assumed to he linear, hut 
time varying. Constraints on hoth the individual decisions 
a,nd on the sum of decisions arc conoid erod. They have 
also derived an algorithm for solving this prohlcm. 

2.5 MUlTICOmODITY FhOh PEOBIjSMS ; 

There are two basic a-pproa.ches which liavo heen employed 
to develop specialized techniques for multicommodity flow 
problems 5 decomposition s-nd partitioning. Decomposition 
approp.ches may be further cha,racterizod as price directive 
or resource directive. 

A price directive decomposition procedure directs 
the coordination between a master program and each of 
several subprograms by changing the objective functions 
(prices of subprograms). The objective is to obtain a 
sot of prices (dual variables) such that the combined 
solution for all subproblems yields an optimum for the 
original problem. There are a nuraber of studies which 
have proposed vafio.tions of this genera,! theme. Tow of 
these are as follows. 

Cromeans, Smith and Tyndall [9] give an approach- 
an extension of the column generation technique used in 
the multicommodity flow problem that simultaneously 
considers network chain selection and ' resource allocation. 
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thus making the problem both manageablo and optimal. The 
flow attained is constra.ined by rosourco a vailability and 
network capacity. 

k resource directive decomposition procedure, when 
applied to a multicommodity network flow pr-oblem having 
K commodities is to distribute the a-rc capacit 3 r among the 
individual commodities in such a. that solving IC sub- 

programs yields an optimal flow for coupled problem. At 
each itero.tion an allocation is made and K single commodity 
flow probleDis are solved. The sum of the capacities 
allocated to an arc over all commodities is equal to the 
arc capacity in the original problem. Hence the combined 
flow from the solutions of the subproblems provides a 
foasihlo flow for the original problem. Optimality is 
tested and the procedure either terminates or a new 
arc capacity allocation is developed. The earliest refe- 
rence of which we are aifsire that suggests, this approo-ch 
appeared in 1956 (Ho backer [4-7]) • However, he did not 
indicate how this approach could he implemented. Sakarovitch[49] , 
Kennington and Shalaby [34] and Held, Wolf, a'nd 
Crowder [25] have presented techniques for implementing 
this approach. 

Partitioning approaches aro specializations of the 
simplex method where the current basis is pe.rtitionod to 
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exploit its special structure. Tiiese techniques are specia- 
lizations of the primal, dual, o]? primal dual simplex method. 
The pa.pers of Saigal [48], Hartman and lasdon [ 24 ] and 
Graves and Mchride [19] deal with primal techniques,’ while 
the work of Grigoriadis and White [2l] is a dual technique. 

ICant [ 32 ] presents an algorithm for finding maximal 
sized sets of flows in a certain class of multicommodity 
flow networks. The class consists of networks with integer 
capacitj?" edges and with each node being a source or oink 
for all but at most one commodity. . This new procedure 
lifts a restriction in T[leitman’s [36] algorithm that 
required flow in edges joining source and sink of the sane 
commoditjT'. 

Golden [13] has developed an algorithm for handling 
nonlinear minimum cost multicommodity flow problems and 
applied it to a. large scale network. The commodities are 
imports and exports j the cost functions considered are 
quadratic and convex. The sotting considered is a Port 
Planning Model which seeks to find optimal simultaneous 
routings through the network while fulfilling requirements 
both at foreign ports and rat domestic hinterlands.. 

The algorithm involves linearizing the cost function 
and solving the resulting linear program, which is, in fact 
a series of shortest route problems. 



2.6 Fixed Cliarge PProblem ; 


The fixed charge problem was first introduced by 
Hirsch and Dantzig [26] in 1954. Since then more tho-n forty 
pa.pers have ajppeared covering various a spec ts of the problem. 
Most of these papers are concerned with the following linear 
fixed charge problem. 

11 


Min Z = + d 


j=l 




s.t. Ax = b (2.5.1) 

X 0 

whore, A is a m x n ma-trix of real scalars 
b is a m vector 
X is a n-vector and 


f 0 
1 

if 

o 

11 

*r 



J 

for all 3 = 1,2,..., n 

if 

x^- > 0 


The quantities d . are the fixed charges which are 
incurred if the corresponding x^ are positive. Eventhough 
more than forty papers have been published in this aroa^ ¥e 
still do not have an exact and efficient computational 
algorithm for linear fixed charge problem. 

In the paragraphs below we summarize some important 
properties of solutions to various fixed charge problems and 
list various solution techniques used, so far. 
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The objective function Z of linear fixed charge 

problem (2.5.1) is concave because c.x. + d .y . is concave 

3 3 3 3 

for ^ 0, and non-negrrbivo linear combiimtion of concave 
functions defined over the same convex set is coucave[22]. 

The constraint set S of (2.5.1) is closed, convex and 
bounded from below. Hence, due to the fact that tho global 
minima of a concove function over a closed, convex set, 
bounded from below, is taken at one or more of tho extreme 
points of the set [ 22]5 tho optimal solution to (2.5.1) 
occurs at one or more of tho extreme points of S. Unfortu- 
nately, for linear fixed charge problem there oxistvS the 
possibility that a. local minima, different from the global 
minima mo.y occur at another oxtreme point. When each 
activity has a positive fixed charge associa-ted with it, 

G "V 0 G la G remo point solution yields a local minima. 

In case problem (2,5.1) has a network structure, an 
extreme point of S corresponds to an acyclic solution, i.o., 
a solution in which arcs carrying flow donot form a cycle, 

Balinski [l] lias proved a number of theorems for 
fixed cost transportation pi'obloms and then on the basis of 
these theorems have presented an approximate procedure for 
the same.. 

Various approaches used sofar for solving linear fixed 
charge problem and its special structures ares 



16 


1. Mixed Integer Programming [44] 

2 . Separable Programming [22] 

3. G-roup Theoretic Approach [45 ]» [46] 

4. Extreme point methods [41], [20], [39], [5] , [lO] , [52] 
[50], [7], [8] 

5. Cutting Plane Methods [5], [5l] 

6. Dynamic Programming [22] 

7. Bianch and Bound Approach [3], [50] 



CHAPTER III 


HETWOSK .ALG-ORITHMS 

In this chapter we will discuss va.rious network 
algorithms which arc used in this dissertation for solving 
subprohloms which a.rise in various solution methodologies . 
"'o will discuss 

1, Mazimal Plow Problem 

2. Shortest Chain Problem, and 

3 * Ilininial Cost Plow Problem. 

3.1 MXBIAL PLO¥ PROBLEM; 

A network can be considered as a pipeline system 
(Hu [30]) with the arcs representing pipelines , the source 
being the inlet of the water, the sink being the outlet of 
water and all other nodes just being junctions between 
pipelines. 

The capacity of each arc can be considered to 
indicate the cross s ectional area of the pipeline with such 
a pipe line system we are interested in the maxiraum flow 
that can be put through it from source to sink. 

In order to make .above intuitive problem precise, 
wo shall first define a flow in the network. A sot of 
non-negative integers X. - is called a, flow in a network if 
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they satisfy the follo'VJing constraints :: 



2 

1 -.- 




i-j- J — s j 

if j Sj t 

if 3 = t 


(3.1.1) 


0;^ 3a (for all i,^) (3.1.2) 

Tho V which appears in (3.1.1) is a non“iieg8.tivo nnjaher 

called the valiic of tho flow, Fote that (3.1.1) represents 

the fact that flow is conserved at every node except the 

source and the sink. Constraint (3.1.2) means that the arc 

flow X. . is always bounded by the capacity of the arc b . . . 

-‘-j 


To find the maxiraTim flow value in any network is 
obviously a linear program with the objective function 
V = 1C 3c_. and the constraints (3.1.1) and (3.1.2). 

7-; s j 

O 

However as this is a very special case of a linear program, 
few algorithms are available which are more efficient than 
the simplex method, which is for general linear programs. 


In the following paragraphs a labelling algorithm 
(Ford and 3?ulkerson [12]) as described by Hu [30] is pre- 
sented . 

The labelling method starts with any arbitrary flow 
(zero flow may be used as the starting flow) and then tries 
to increase the flow value. The algorithm terminates if the 
flow value can not be increased. The algorithm is a systematic 
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way of searcliiiig all possillo flow augraenting paths from 
tig to 1.1.. This is done hy giving labels to nodes to 
indicate the direction that flow in an arc zaay be increased 
once a flow augmenting path is found. The flow is increased 
to its mariraum capa.city a,long the path and all the labels 
on the nodes a.re omsed. ¥e then start to assign now 
labels to nodes based on the now flow. 

The algorithm consists of two .steps z Step 1 is the 
labelling process, which assigns labels to nodes, and 
Step 2 is flow change. Steps 1 and 2 are iterated until 
the incre8.se of flow becomes impossible. 

STEP 1; 

Labelling Process. Every node is always in one of 
the throe states, la.bollod and scanned, labelled and 
unscanned, or unlabelled. A node is labelled and scaained 
if it has a label and we h8.ve inspected all its neighbouring 
nodes (Two nodes are neighbors if both are connected by 
one arc. ) A node is labelled and uiis canned if it has a 
label but not all its neighbors have been inspected. A node 
is unlabelled if it has no label. 

Initially all nodes are unlabelled. A label for 
node IT- always has two parts. The first part is the index 

tJ 

of node which indicates that xfl-e can send flow from 

N. to E • , and the second part is a number which indicates 
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the maxii'ium amount of flow wc can send from source to IT . 

3 

without violating the capacity restrictions. ¥e first 
assign the lahol [s"^, 6 (s). = CiO ] to 11^. The first label 
simply says that wc can send flow from to itself | the 

uJ 

iiumbor indicatos that there is no upper bound on how 
much can be sent. The node is now labelled and unscanned 

o 

and all other nodes arc unlabelled. In general, select a 
node M - which is label3.cd and naunscannod , Assume F- has a 

w J 

labo3- of the form [i"^, ^ (j)] oi" [i"”!- To all 

neighbouring nodes of IT. which are unlabelled and for 
which O^x-i^^b.-,.. assign the label ^ (k)] to 

where 

€. (k) = min [£( 3 ), b^|^ - (3.1.3) 

The + and - signs in the first label indicate how 
the flow should bo changed in Step 2. Fow all tho neigh- 
bouring nodes of F- have labels; F- is considered to be 

J J 

labelled Sund scanned and may be disregarded during the rest 
of this stop. (If one inspects all the neighbors of IT- and 
can not leibel the neighbours, theai a3-So H- is considered to 

J 

bo a labelled and scanned node). All the nodes I'Tj. a,ro now 
labelled and unscanned. 

Continue to assign labels to hneighbours of labelled 
and unscanned nodes until either is labelled or no more 
labels ca,n be assigned and is unlabelled. If F^ cannot bo 
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labelled, no flow augmenting path emists and, hence, the 
flow is ma.simal. If-U^ is la,belled a flow augmenting path 
can bo found using Step 2. 

STEP 2; 

‘Plow Changes Assiuac that llj. is lEibellod [K"^, ^(t)]. 
Replace by + £ (t) and turn to If is 

labollod [f,€{k)], replace by Xj, H-£(t) and tune 

to N-. If ITt^ is labelled [j , S (k)] rcialace x, . by x, . ~6.(t) 

J -K- xj xj 

and turn to If_- . Continue until H_ is reached. Erase the labels 
J o 

on all the nodes and go back to Step 1. Pew other algorithms 
are also available for the same problem, such as Dinic [ll], 
ICarzanov [33] and Malhotra, Kumar and Malieshwari [38]. 

These algorithms have better computational bound tlmii the 
original Ford o,nd Fulkerson's [12] algorithm. 

Label].ing algoring described above (Ford and 
Fulkerson [12]) has been coded into FORTRIF IV for the 
computer system IBM 704v/l4Ql. The listing for above code 
can be found in Appendix B. 

3.2 SHORTEST GI-IAIl PRO BLM; 

In this section we shall discuss one of the most used 
algorithms in this dissertation for finding shortest chain 
between any pair of nodes in a network. The problem of 
finding shortest chains arises very frequently as a sub- 
problem in other problems of optimizati on and also has a 



practical interest of its oim. Every arc of the 
network has associated with it a dista.ncc The 

problem is to find a chain from any node to any other 
node N-, with the sum of the distances c. . of o^rcs in the 
chain a minimum. 

The algorithm used b 3 r us for finding shortest chain 
between any pair of nodes in a network is called as 'revised 
cascade method' (Hu [28]). Tho matrix C is a iT x IT matrix 
(H is number of nodes in tho network) onid is called as 
distance matrix. At the termination of a,lgorithm c. - 

X J 

represents length of the shortest chain from node i to node j 
C matrix provides lengths of the shortest chains at tlie 
termination where as route matrix E provides the routes for 
shortest chains. 

The algorithm proceeds as follows. Initially all 
c. - entries are direct distances between node i and node 3 ./ 

J- 

c^j^ « 0 for i~lj... 5 T'T. 

c. . = ^ If no directed a.rc joining IT. to IT- 

a 3 ^ w 

is present. 

The initial R matrix has entries of r^j^ = k for each 
nonzero entry. 

How let j bo fixed at j = 1,2,..., IT. For each value 
of 3 , say 3 = 3 ^ repeat the operation 3 . 2,1 for each entry 

°ik i 5^ Oq 
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am 


'ik 


ij 


+ c 


ik] 


(3.2.1) 


TliG matrix thus geiiGrated is shortest distance matrix. To 
deteriaine tiie shortest routes, a route matrix is generated 
while the shortest distance matrix is being constructed. 



ifc., •<" c-. + c., 
ik ^ ij jk 


r- ■ otherwise 

-L. J 


( 4 . 2 . 2 ) 


This completes the description of revised cascade 
method for finding shortest chains. This raethod was coded 
in lORTEAlf IV for IBM 7044/1401 system. The listing for the 
same can be found in Appendix B. 

3.3 MIKIMI^ COST PM PEOBLEIk 


Hu [2S] has given an algorithm for solving minimura 


cost floxf problems where the shipping cost over an arc is 
a, convex function of the number of units shipped along 
that arc. The problem can be sts^rted as 


Min Z = I 


Id 




(3.3.1) 





for d = s 

for j / s,t (3.3.2) 
for d = 


where c . ^ ) are nonnegs-tive convex functions of x . • , 

(c. .(0) = 0) and the arc flows x. . are required to bo 
Id ^d 



which is 
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positive integers or zero, Tho paraiietcr v, which is 
required to he a noimogativo integer, represents the total 
flow from source to sink. Equations (3.3.2) express the 
conservation of flow at nodes other than soured and sink. 


AhQORLTmi: 

The a3.goritha for solving the 'aininua cost flow 
prohlem in convex cost networks can he sinplj'' described as 
follows. Starting with all x. . =0, send one unit of flow 

from Eg to along the path wdioso incroncntal cost relative 
to the existing flow is ninimun. 'This can be done "oy any 
of the existing shortest-path methods with up costs and 
down costs as lengths.- (Up cost and Dowm cost for a directed 
arc are the cost of sending one additional unit in the direction 
of the arc and opposite to it respectively). In the ^ 

beginning only up costs are relevant since there are no 
positive arc flows which could be potentially cancelled. 
Redefine the up costs a,nd down- costs based on the now flow 
pattern obtained, and send one additionad. unit of flot-r along 
the path -with niininui-a incremental cost . The process of 
using the minimum incremental cost path is repoatod until 
the total outflow of is v (or the total inflow of is v). 

This algoritlim in principle is used in the heuristic 
methods developed in Chapter IV. This algorithm is also 
coded in P0RTEA.N IV for IBM V044/1401 system and is 
presented in Appendix B. 



CHAPTER IV 


CAPACITY EiCPANSIOH IN SIIGLE COffl-IODITY FLOW fJBTWOEES 

In this chapter wc will discuss -various heuristic 
methods and an exact method developed for single commodity , 
flow netifork capacit3'' expansion problem^ before that we 
will discuss various types of capa.city expansion costs. 

4,1 CAPACITY EXPAHSIOH COSTS' 

Capacity of a network cam be increased in two ways, 
either by adding some arcs to the network or by improved 
maintenance of arcs. This may well be explained by con- 
sidei-ing a railway network. The former case is q,uite 
obvious because addition of new links x-Jith some capacity 
will naturcally improve ca.pacity of the netxrork. In the later 
case improved maintenance of the links in existing network 
will enhance trains to operate at a faster speed resulting 
in increased capacity of the rail network. 

Capacity expansion cost function will naturally 
depend on the way in which capacity expansion has been 
achieved. Some of the cost functions which arise in prac- 
tice are discussed beloxf. Fulkerson [16] has assumed the 
cost of increasing capacity of a link as linear a.nd 
homogeneous. 
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If capacity of the network is increa,sed by adding 

3-inks then total cost may look like something in Fig. 4.1.1. 

The sudden jumps at some interval of . represent the cost 

-^3 

of building a new link where o-s vs-riation in cost in o-ii 
interval of x. is because of imintenance and shipping costs. 

-t, J 

If coi-pacity of the network is increased by adding 
new links and if mainteimnce costs can almost be neglected 
then the total cost function may look like something in 
Fig.- 4.1.2 or Fig. 4.1.3. In. Fig. 4.1.2s shipping costs are 
assumed to be concave where as in Fig. 4.1.3 shipping costs 
are assumed to be linear. 

If capacitjr of network is increased by adding new 
links and if shipping costs are negligible as comioared to 
costs of building new links anad maintenance cost then the 
total cost function may look like something in Fig. 4.1.3 
of Fig. 4.1.4. In F ig. 4.1.3 it has been assumed that . 
maintenance costs are linear where o.s .Fig. 4.1.4 incorporates 
more common type of maintenance cost functions the convex 
function. If capacity of network is increased by improved 
maintenance of links only then total cost function will u 

look like in Fig. 4.1.5. 

In this disserta-tion we have assumed that maintenance 
and shipping costs are negligible as compared to costs of 
building new links . Thus total cost function whicli we 
refer in this work will look like in Fig. 4.1.6. 
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4.2 MTHEMA.MCAL FOBIOTMTIOFi 


Noinonclature ; 


% 

^±2 


1 

b. . 

d?r . 

10 


^iO 

Y 


1 

^iO 


Mode identification number 

Directed arc joining and and heading 

towards M. 

Arc ca,pa,city after 1-th levoO. cost is incurred 

l_i 

Cost of increasing ca,pacity from b. to 

J 

Plow on A. - 

-L J 

E.equired flow from to 


0-1 Ya.riable. 


1 

^iO 


1 if x_. ^ ^ b^"^ 


10 


0 Otherwise 


■0 


can be i-epresented mathematically as follows. 

Object i ve Punctio n; 

The objective function is to minimize the total cost 
of introducing va.rious arcs o.t vo.rious levels in to the 
network. Cost of introducing Ic-th level of can bo 

expressed as d^yy^-. y^- is 0-1 variable. It should be 


noted that A| ■ ca-n be introduced at k-th level onl 3 r and' 

only if A. . exists in the network with k-1 levels. The 
10 

total cost function, TO, can be expressed as 


TO 


^ ^ ^ 

i 0 k 


,k k 
d . . y . . 
10 ^10 


(4.2.1) 
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The objective function v.2.1 is to be optimized 
considering the following constraints. 


Constra ints 2 


Thero are two type of constraints for this problem. 
The first type of constraints ai’O Xirchoff's Hodal Equili- 
briuTii constraints wherea.s second tsrpe of constrs-ints arise 
to enforce the y^ . variables 8 .s 

1 tJ 



1 


ix 





0 Otherwise 


( 4 . 2 . 2 ) 


Kirchoff^s Fodal Eg ^uil ibri u m Cons t raint s ; 

These constraints essentially represent the fact that 
flow is conserved at every node except the source and the 


sink. This can be represented mathematically for node J as 


/ -V if j = s 

X X. . - X X., = J 0 if 3 s, t (4.2.3) 

i k \ 

V If j = t 

These constraints will be E in number. 

It should be noted that convent! ona3- capacitor cons- 
traint of minimal cost flow problem is implicit in this ^ 
formula. tion with some variation. In minimal cost flow 
problem arc ca,pacity can not be increased i.o. cost asso- 
ciated with increasing the ca,pacity is infinite, where as 
in this formulation arc capacities ca.n be increased by 
another level by incurring some finite cost. 
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OtliGr Constraints 


These constraints enforce variables a-s per 4.2.2. 
These constraints can be expressed mathematically as, 


1—1 1 
X. . - b. < M j. . 
13 


(4.2.4) 


This type of constraints xirill be as many as there are y^^ 
varia-bles. 


4.3 BmiTCH MID BOUND I-IETHOD (EXACT 14BTH0D),: 


The problem (4^2.1) - (4.2.4) can be solved, by using 
any of the existing 0-1 programming algorithms. 

The exact method which we are a.bout to discuss in 
this section is essentially a tree search method. Bor this 
method there need be no specific relationship between level 
costs and level capacities for all arcs. 

The method starts with an upper bound set by input 
feasible solution obtained by any one of the heuristic 
methods of Section 4.4. If a feasible solution is not 
available- an arbitrarily large number is sot as initial 
upper bound, however a good inpiit feasible solution may 
result in sizeable saving in computational effort. 

The branch and bound approach is just the simple 
*divide-and -conquer’ strategy. Roughly speaking it repeatedly 
partitions the sot of all feasible solutions in to smaller 
and smaller subsets, termed 'the branching (partitioning) 
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process’, and for each of the subsets, a lower ’bound for 
the cost of solutions within the subset is computed, termed 
the 'bounding process'. All the subsets are listed in a 
list called the 'master list'. The best known feasible 

solution is the 'incumbent' solution. After each parti- 
tioning a subset for which, either the associated lower 
bound is the lowest lower bound, or it ha,s been ascerta-ined 
by some means, (e.g., if it is empty, or if the associated 
lower bound is greater than 03.- equal to the objective 
function value for the incumbent solution), that it can not 
contain a feasible solution better than the incumbent 
solution, is discarded from all further considerations, 
(i.e. , deleted from 1 ) termed a.s the subset having been 
'fathomed'. The process is continued in this fa,shion till 
either the master list ' 1 ' is empty, or the loirest lower 
bound is close enough to the objective function value for 
the incumbent (incumbent value). On termination incumbent 
is an optimal solution. 

In the branch and bound method discussed below 
branching takes place whenever lower bound on total cost 
is greater than or equal to the upper bound . 

Whenever a feasible solution vrith total cost less 
than the present upper bound is found it sots up new upper 
bound on the total cost. 



The following is a stop wiso description of the 
algorithm. 

1. Find a teo.si'ble solution to the prohlon using a.’iy 

of the houristic methods discussed in next soction. The total 
cost associa,tod with this solution will sot oni upper hound 
IJ for the prohloj'i. 

2. liumhGjr the arcs (bi = 1,2,... , M). Eoxli arc cooa he 
added to network n,t one of t possible levels , 

3. Let a,ll arcs ho present in network, each one at its 
first level. 

4- Sum up all individual arc costs for present configu- 

ration to co.lculato total cost D. In the process of summing 
if at any timo sum exceeds upper hound go to Stop 7. Otherwise 
after ca.lculo.ting D go to Stop 5- 

5. Find out maximal flow g for present configuration. 

Go to Stop 6, 

6. If g is gi'eater thorn or equal to required flow v, 
replace previous solution which set upper hound hy present 
solution a.nd go to Step 7. Otherwise go to Step 14. 

7. Set a = 1. 

8. If M - a is less than or equa,l to zero go to Step S, 
Otherwise go to Stop 10. 

9. Print prosont stored solution. This is optimal 


solution. Stop. 
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10. Sot IncremGnt level of h.--th arc liy 1. let 

new level for li-tli arc be f. Go to Step 11. 

11. If f is greo-ter than t : go to Step 12 ^ otherwise 
go to Step 13. 

12. a -^—a+l, go to Step 8. 

13. Set level 1 for arcs h+l to H (both inclusive). Go to 
Step 4. 

14. Increment level of M-th a.rc by 1. Let now level for 
H-th arc be f. Sot h-^ II, Go to Step 11. 

Il lus t r ation s 

Consider following problera. 

Fig, 4.1,2 


Number of levels considered for each arc = t =2. 


Arc Arc 

? I 

1st level 
cost 

1st level 
capa,city 

2nd level 
cost 

2nd level 
capo, city 
(TotoJ. co.p.ac 

1, 2 

1 

o 

CO 

5 

3.0 

10 

1, 3 

2 

2.0 

7 

7.0 

12 

2, 3 

3 

5.0 

5 

2.0 

11 

3, 4 

4; 

7.0 

5: 

o 

CD 

17 







Pig. 4.3.1 Plow Digram for tHo Branclx and Bomid Algorithm 
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Let the required flow fron node 1 to node 4 he 10. 

Step 1 Since we do not have a feasible solution with us lot 
initial upper bound on total cost be 100.0. 

Step 2 Arcs are numbered as shoTO in the figure. ■ 

Step 3. Let all arcs be present in the network o-t their 

first level. Low we arc at node A in the figure 4.5. 

Step 4. Total cost D = 22. Since it is less than upper bound, 
we go to Step 5. 

Step 5 Hazimal flow g = 5. 

Step 6 Since nazimal flow is less than required flow this 
is not 8. foo,siblo configura.tion. Hence go .to 
Stop 14 . 

Stop 14 Increment level of 4th, arc to f = 2. Low we arc at 

node B in the tree of figure 4.5.3. He go to Step 4. 

Step 4 Total cost D for this configuration =50. 

Step 5 Maximal flow g = 10. 

Step 6 Since g = v new uppei' bound on total cost = 50 •= U. 

Step 7 a = 1 

Stop 8 Since li - a = 4-1 = 5 is greater than zero we 
go to Step 10. 

Stop 10 h = 5. Low level of 5rd arc = f = 2. 

Stop 11 Since f ^ t go to Step 15 . 

Step 15 Set level 1 for arc 4. How wo are at node c of tree 

of Figure 4.3.3- 
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Step 4 ‘Total cost D = 24. Sinco D<C U we go to Stop 54 

Step 5 Masiao.! flow = g = 12. 

r 

Step 6 Since g v new upper bouiic'. on total coot = 24 = U, 

VJe proceed in this way until all the nodes from A to 

P of the tree in Figure 4*3.3 are iniplicitel 3 :" enumerated. 

4.4 HEURISTIC METHODS; 

Single commoditj?' capacity Gxpa.nsioii (SCCEP) problem 
ca,n be .solved by using the exact algorithm of Section 4,3 or 
possibly by some 0 “ 1 programming algorithm. But for the 
exact method of Section 4,3 computations required increase 
at an exponential rate with the number of arcs in the network, 
hence it may become time consuming for networks with large 
number of arcs. Honco an efficient algorithm for solving 
the Single Commodit 3 '‘ Capac it Expansion Problem is needed. 

Here we present four approximate solution techniques for 
solving single commodity capacity expansion problem which 
in the absence of a.n efficient exact algorithm may be proved 
useful. 

Nomencl at ur e ; 

I'Tj^ Node with index i 

Source node 

O . ■ 

Ki. Sink node . 

A - ^ Arc joining N.; and N- and directed towards N . 

10 ^ 1 j J 
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d - . 

13 


Cost of introducing JL-tli level of in 

the notuork. 


P' 


10 


Jr 

d : . 

10 


k 

^iO 


“^iO 


capacity attained by A. • o.ftor having 

10 

incurred d^j 

Cost of introducing final level (k-th lovel) 
of A- - into the network. 

ti 

Co.pacity attained by A. ■ after having 

-1- tj 

incurred d^j 
Capacity of A 


iO 


(b. - will nonnall3’' be equal to p,- ^ ? 1 = l,,..k) 


Heuris t i c__-_ 0. 


The a.lgorithra first finds out single minimal cost 
chain from iTs to which ha,s got capacity to hold a flow 

groa.tor tlmn or equal to required floa' v. It should be noted 
that the algorithm starts with an self imposed constraint of 
finding single drain which can meet flow requirement. In 


the lator stops of 8.1gorithm arc in the chain having largest 
final level cost is co.pr'.citatod to next lower level and 
cost per unit flow is upda.ted a.ceordingly. Other arcs in the 
chain are capacitated to the capacity of the existing fina,l 
level and cost per unit flow is set to zero. This procedure 
is repeated for all arcs in the chain with their fiiml level 
costs in descending order. 



40 


The above mentioned procedure for capacitating arcs 
and modifying c-^'s makes an attempt to remove the self 
imposed constraint of he.ving the solution as single chehn. 

The following is a stop wise description of the algorithm. 

1. let required flow from source to sink bo v. 

2. Assume that on each A- • there is flow equal to V. 

X J 

3. Calculate cost per’ unit flow c . - for each A. . as 

X J X J 

follows. 

c. ■ = Cost incurred to aquire arc capacity greater 

tha,n or equal to v/v. . 

4. Find tlio shortest path between and the with 

as distanoo between If- and N-. let F denote the set 

tj 

of A-^ constituting the shortest path, (let total no. 

of A. •£. F be u). This gives first feasible solution 
, ^ 3 

as a single chain. The total cost corresponding to this 
solution sets up the upper bound U on total cost (U = l). 


5 • Set q = 1 , 

6 . Sort A. . C such that dA^. are in descending order. 

XiJ 

7, Set w = 0. 

Calculate new and impose now for all 

as follows, 

a. If A. . F and ranks q-th in the ranking done above 

-K ^-1 

- . : : then, hij = : 

k-1 : T 

^ . dt 


1=1 

c . . = — 


13 
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^13 


1r 


P 


13 


13 

k~l 

-i3 


b. 


1' ancl does not rank q~th in tlio ranking 
done above. 


13 


'^13 


“ij = 0 

If . ^ F 


k 


c . . = 

3-3 


1=1 


d^ . 

13 


V 


8. Using Hu’s approach as doscrihed in Section 3.3? find 

the cost of aquiring the required flow with c. • as 
the cost per unit flow on with capacity restric- 
tion K-.-. Also find flows on each A. . this is one 

X J X J 

of the solutions. If tot.al cost D' corresponding to 
this solution less than U store this new solution 


and Set U = D' . 

9. Incromont q hy 1 (q<^q+l) ■ 

10. If q is greater than u, print out storod solution and 
the value of U which will he total cost value for this 
solution. 

If q-^ u go to Step 6. 

The flow diagram for this method is as shown in 
Figure 4.4.1. . : ' 
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H_ouristi c _ 2 ; 

Heuristic 1 ntukes an attempt at each, iteration to 
removo one level of an arc which gives maximuin decrease to 
the objective functiono Once a level of an arc giving 
maximum decrease in objective function is removed ?ao further 
attempt is made to improve this solution. Hoxt iteration 
sta.rts with next arc in the very first solution giving 
better docrease in objective function. 

Heuristic “ 2 differs with Heuristic 1 in the sense 
that it makes some attempts to improve the solution obtained 
by removing ono of the levels of an arc in the very first 
solution. 

The following is a stepwise description of the 
algorithm. 

First 7 stops for this method are same as in 
Heuristic 1, 

8. Using Hu's approach as doscribed in Section 3.3 find 

the cost of acquiring the roquirod flow with c. . as 
tho cost por unit flow on il- - 1 x 11-1 csipacity rootric- 
tion h.v. Also find flows on oo-cli This is ono 

of the solutions, let G be tho set of A. . having non 
zero flow on them. If DV- the totrul cost corresponding 
to this solution is less than U set U_= 1* and store 


now solution 



43 


Find A., sucli thnt A- - 6 G and dt is largest. 

1 J -L J J 

Calculate new c- and impose , neu b. -’s for aJ-l 

X J X J 

A. ^’s as follous. 

X J 

s 

a. If ^ G and eq.uo-ls A^^^ 


then 


k -1 

Pi1 


ij 


k -1 


I 

b. If A. 4X G and does not equal k. ■ 


T f A • - G 

J.X 


^ dv, 

ill 


^ If w^4 go to Step 12. OthorwiSG go to 

Step 8. ■ ■ 
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12, . Increment q ly 1 (q^r*’ 1+1 ) 

If q ^ u print tlie v.alue of U, wliicb; will be tlie 
tote.l cost corresponding to the best solution 
ohtained, and the stored solution. If q u, 
go to Step 7. 

The flow dia.gram for this method is ns Bliowj;i in 


3J'igure 4.4.2. 


Heuristic 3 « 


This heuristic method is almost same ns heuristic 

*s 

method 1 except the criterion for sorting ^^ch 

tha,t 5 


j 


d j 




n" 


.k~l 


It was felt that- capacita-ting an arc in the first chain 
solution to its next lower level on the basis of large 
will some tines be iiioro appropriate than raeroly largo 
following illustration will make this point more clear. Let 
two arcs in the first chain solution ha-ve the cost structure 
as in Figure 4.4.3. There is flow of 18 unite of each arc. 
Arc 1 has got its final level cost as 12 (23-11) c-nd its 
final level has a flow of 10 units (18 - 8).: Though arc-2 
has its final level cost as 10 (28 - 18), which is lower 
than that of arc-1, it carries only one unit of flow. 
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Arc— 1 Arc— 2 

Pig. 4.4.5 


It mGaiis for arc - 2 additional cost of 10 is incurred only 
for one unit of flow^ where as arc - 1 incurs extra cost of 
12 for 10 units of flow. Hence, it will seem more logical to 
capa.citate arc - 2 at its next, lower level than arc - 1 xHiicIa 
otherwise will be capacited to its next lower level according 
to criterion of heuristic - 1 




"^■'5 Gompu ta ti onal Perfonnance ? 

The algorithms described above have been coded and 
tested on a number of randomly generated problems. Three 
problem sets were solved by each of the heuristics. The 
computer programs were written in TORTPufi-F IV for the 
IBII 7044/1401 system and were run using the P0RTRA.N compiler 
of that machine. 

All test problems were generated as follows; 
starting with the required set of n nodes and an incidence 
matrix for the network , level costs and level capacities 
for each arc were generated randomly from a uniform distri- 
hution. Total number of levels generated for each arc were 4 
For each node size three test problems were generated. 

Table 4.5.1 shows the best solutions given by various 
heuristics for problem set no. 1. It cam he very easil3/' 
observed tha.t heuristic-2 a>,lwa3rg gives a solution at least 
as good as that given hi?' heuristic -1. Sa/ne holds true for 
heuristic-4 and heuris tic-3 respoctivel.v. 

In many of the cases tho best solutions given h:y 
heuristic-2 and heuristic -4 are not mueli better than that 

of houristic-1 and heuristic-3 respectively. Considering 

/ 

the time taken bj?" heuristic-2 or heuristic— 4 to solve the 
large problems it seems to be a good deal to solve large 
problems either by heuriStic-1 or houriacic-3. 



But if nearness to the optimal solution even by a. 
small amount is important heuristics 2 or 4 seen to be the 
desira-ble methodologies. Table 4.5.2 shows the average 
time taken to solve 3 sets of loroblens of various node 
sizes, and average number of iterations for getting the 
best solution given by the heuristics. Figure 4,5.1 shows 
the plot of average tine for three sets of problems takeii 
by heuristics 1 and 3 V"S. number of nodes in the problem 
on linear scale. Figure 4.5.2 shows the plot of average time 
taken for three sets of problems 'bj heuristics 2 and 4 vs, 
number of nodes in the problem on linear scale. 

Bjr observing the curves of Fig. 4.5.1 and 4.5.2 it 
wo.s felt that some polynomial can be fitted in to these 
curves hence these data, for heuristics 1 and 2 were plotted 
on log- -log paper in 3?igure 4 . 5 . 3 . For both the heuristics, 
the plots can be capproximated by straight lines. These straight 
lines can be represented as log T = m log .4 -!- log C oi’ siniply 
T = where m is the slope of the straight line in Fig. 4.5.3, 
C is a constant, F is number of piodes and T i.j computatioiia.l 
time. For heuristic-1 coiaputational tine can be i-oprcsented 
8.S y 

T = 7.6 X 10“^ x 

For huuri&tic-2 computational time can be represented as 
= 7.45 X 10“^ X 
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Because of erratic functioning of the computer system 
it was not possible to solve many problems with bra.nch and 
bound method. However, for some of the problems solved 
computational time and optimal solution is shown in 
Table 4.5.3 along with computatioiial time te.ken and best 
solutions obtained by ve^rious heuristics. 

Percent deviation from, the optimal solution for the 
solutions obtained by various heuristics are plotted in 
Pig. 4.5.4 and Pig. 4.5.5 against node size (both on 
linear scale). 

Prom Pig. 4.5.4 and Pig. 4.5.5 it can be observed 
that heuristics - 2 and 4 provide solutions quite near to 
optimal solution than heuristics 1 and 5 respectively. 

Pig. 4.5.6 shows the time taken by branch and bound 
method to give optimal solution for various node size 
problems. 

Conclusions and Scope for Further Research; 

It is quite difficult to solve the problem considered, 
optimally. The computational time rises steeply vrith problem 
size. Even for the linear fixed charge problem which can be 
thought of as our problem simplified, the available exact 
methods are very inefficient. 
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For tlie type of problem considered very limited work 
is available in the literature. Following avenues for further 
research are suggested. 

1. Efforts need to be made to solve the problem consi- 
dered efficiently. No efficient and exact method is ave.ilable 
for the problem. 

2. Multicommodity flow problem with cost functions 
considered in Fig. (4.1.1) - (4.1.5) has got resemblance 
with many real life problems and almost nothing has been 
done in this area. Thus efforts can be made to solve 
multicomiaodity flow problem with these cost functions. 
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Table 4.5.1 


Number 

of 

Nodes 

Best solutions given by 



HcuristiC“l 

Heuristic “2 

Heuristic~3 

Houristic--4 

3 

306.97 

306.97 

306.97 

306.97 

4 

364.79 

350.68 

364.79 

361.94 

5 

508.12 

485.57 

508.12 

458.39 

6 

358.46 

253.49 

259 .54 

256.43 

7 

416.5 

390.78 

438,18 

342.60 

8 

534.22 

492.59 

465.64 

465.91 

9 

308,41 

285.41 

526.47 

300.41 

10 

63S.32 

318.25 

318.25 

310.4 

12 

777 . 68 

656.7 

904.32 

697.4 

14 

436.39 

427.50 

489.23 

415.34 

16 

603.71 

592.55 

663 .44 

490,54 

18 

618.39 

551.49 

689.77 

586.26 


IJ.T. KrMKTUR 

t. la RARf 

: ' “STM?. 
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Tabic 4.5.2 


ITumbor 
of nodes 
in net- 
work 

Hcuri 

,stic~l 

Houristic-2 

Heuristic-3 

Heuristic-4 

Time"*' 

Itera- 

tion 

Tirae'*' 

Itera,- 

tion 

Timo'^ Itera- 
tion 

Time'^ 

Itera- 

tion 

5 

0.5 

2 

0.5 

2 

0.5 

2 

0.5 

2 

4 

1.5 

2,5 

1 

3 

1 

■2 

1.33 

2.67 

5 

1.33 

2 

2 

2.5 

1 

2.5 

1.33 

3 

6 

1 

2.33 

2.5 

3 

1 

2,5 

2.33 

3.33 

7 

2.33 

2,33 

6.3 

2.67 

3 

2.5 

6,67 

7.3 

8 

3.25 

2.25 

■ 10 

3 

3.5 

2 

9.5 

4,2 

9 

5 

2 

17 

5 

5 

2.75 

14 

3 

10 

8 

4 

21.5 

9 

6.5 

3 

21.5 

6.5 

12 

16 

2 

37 

6 

15 

2 

35 

4 

14 

19 

3 

63.5 

7 

19 

2 

64 

6 

16 

28 

4 

101 

11 

28 

4 

100 

7 

18 

36 

5 

192 

5 

37 

3 

192 

14 ^ 


+ Average time is in IBM 7044/1401 seconds for PORTRiN 
compiler. 



53 






LA 

VO 













CO 

O 







d 







• 







0 


' 1 

D^. 

o 

KV 

LA 


DA 






•j— 1 


^ 1 



iH 

H 








-P 


* 












d 


■ o 












J — [ 


•H 


i>- 


CA 

DA 

O 






0 


n 1 

rh ® 

fS3 • 

. (j\ 

• 

CA 

n 

DA 

« 


MO 

'ft 






0 


' -H i 


VO 

H 

00 

VO 

CM 






H 




o 

VO 

LA 

LA 







cd 


! 3 



DA 


CM 

DA 






a 


? -o ! 












*H 


iii 1 

EH . 

LTV 

• 

H 

rH 

CM 

VO 






-P 

H 


■ ■■ i". 


O 










0 


i 



Ai 











, 1 



CO 


CM 

t-- 






a 


1 



o 

« 

DA 

ft 






0 


j 

CJ^ 

o 

DA 

GO 

O 

H 






H 


1 

1 



H 

CM 

LA 

DA 








KM 



CA 

CM 


00 






0 


II 




r-| 

LA 

r-i 






•H 


Of 


0 

« 

* 

O 







M -3 


•H| 


VO 


00 

(A 

00 






0 




o 

VO 

o 

LA 

CA 






•H 


CQI 


{K\ 

DA 

LA 

CM 







H 


»H I 












d 




LTV 

rH 

rH 

H 

DA 






0 




• 










rd 



EH 

o 




00 








Ixli 



UD 

[>> 

CM 







|>v 

1 





• 

00 

« 






40 





» 

CM 

• 

JA- 




'"<J 




CMf 

' 1 

o 

a'^ 

CM 

CM 

Hi 




O 


•TO 

1 

1 1 










M 


0 


Ol 


t> 

CO 

o- 

OA 

CD 




40 


d 

1 

•H 1 

J 

cn 

VO 

LA 


C-- 




0 


'H 

1 

H^! 


® 

« 

'« 

0 

• 




a 


cd 

i 

ml 

tS 3 ■ 1 

VD 

o 

LA 

DA 

o 






-}0 

1 

‘HI 


O 

LA 

00 

LA 

CA 




'TO 


rCi 

s 

HI 

% 

IK\ 

A 

'M- 

CM 

CA 






0 

i 



LA 

rH 

CM 

CA 

^ , 




0 



^ 1 

Ml 

eh 

e 




\ 

j 




40 



• 




Ai 



1 






d 





00 


{> 





'TO 

0 

0 

i 

1 


i 

o 


00 

* 




d 

•H 

•H 


1 

CJl 

r O 

DA 

GO 

• 

LA / ' 




cd 

40 


»! 

HI 


j 

jH 

CM 


CM 





0 

d 

O j 

3 1 










M 

*H 

rH 

H 

OJ 


! 1>- 

0'^ 

CM 

VO 





0 

H 

0 

rQ , 

>H I 


O*^ 

t- 

rH 


LA 

M 



d 

3::3 

0 

CS If 

-P 3 

^xi : 

0 

• 

o 

0 

» 

H 



cd 

0 


Eh " i 

cal 


VO 


00 

CO 

VO 

O 

H 


H 

""M 

43 

i 

*H I 


o 

LO 

o 

LA 

H . 

' 

o 




0 

i 

.HI 


CA 

DA 

LA 

DA 



jS: 




0 

1 

^::il 







a> 

+:> 


>5 

>> 

40 


Ol 

M! 


LA 

<0 






0 

d 


r-Q 



1 



O 

rH 

H 

rH 

CM ; 

CD 



d 

•TO 

0 

1 








H 

0 


0 

0 


1 








-P 

rd 


f> 

d 

d 

i 









•P 


'H 

•H 

0 

1 




O 

iH 

DA 

O 



© 


cd 

•H 

j 



CA 


tr- 


VO 

•H 

d 

El 


-p 

-P 

1 

. O 

• 

« 

o 

0 

e 


•H 

•H 

d 

rQ 

cd 

! 



VD 

O 

LA 

VO 

CM 

CQ 


-P 

0 

0 

•H 




O 

CM 

0"v 


CA 

0 

m 


•H 


> 

1 



^A 

DA 

DA 

Csl 

DA 


o 

H 

-P 

d 

0 

i 








O 

H 

cd 

d 

0 


i 








d 

cd 

d 

rH 

‘H 


. } 










o 

0 

-P 

CD 

\ 

o 









-H 

ca 

d 

fclQ 

I 

n 







d 

o 

-p 


rH 

Cd 

i 

03 

'TJ 








cd 

H 

0 

-p 

! 

H 







H 

H 

-p 

cd 

IQ 

d 

■ f' 

rQ 


AJ 

A 

CM 

VO 

DA 

0 

0 

d 

a 


0 



H O 



CM 

CM 

am- 

rQ 


H 


4 D 

0 


, ■ 

cd M 





CA 

§ 


a 

-P 

0 

H 











o 

H 

0 

0 










M 

o 

0 

M 



' Hi 

1 

K^ 

LA 


00 

O 









1 





H ■ ■■■.:■: 

11 

II 

11 

11 

11 

11 

i 


:-i 

DA 


LA 

VO 

■CA 

d 

H 

EH 

0 

bO 

i:<3 







^ HEURISTIC-2 
HEURISTIC- 4 


NUMBER OF NODES N 

2 PERFORMANCE OF HEURISTIC- 2 AND 
HEURI5TIC-4 


















inj' AL TIME IN SECONDSCLog 












u ‘S ',"" 




• HEURISTIC-1 
^ HEURISTlG-2 


2:' "4 6 8 12 16 20 

NUMBER OF JOBS N 


NUMBER OF NODES FOR HE 
AND HEURlSTIC-2 






^■s 




■ 










COMPUTATIONAL TIME IN SECONDS 








200 f- 


100 












: NUMBER OF NODES N 

FIG.4-5-6 PERFORMANCE OF BRANCH A 




METHOD 







CHA.PTSR ¥ 


CAPACITY EXPANSION IN I^UlTICOffilOPITY PLOW NETWORKS 


5.1 INTRODUCTION; 


The Maximal flow prolDlem involves finding maximum 
flow from source N^ to the sink N^ without violating the 
capacity constraints x^^^ ^ on the arcs. If there are 
many sources and man37' sinks and if the flow from any source 
can be sent to any sink, then this problem can be converted 
trivially into a one source and one sink problem by creating 
a super source and a supersink. If we make the restriction 
that the flow from certain sources must be sent to certa.in 
sinks, then the problem is a, multicommodity flow problem. 

Lot there bo sources IK and sinks IT , (s = 1,..., q j 

O o. 

s' =l',...,q'), where the flow s is from Ng to N^ , . Let 
x?^ be the flow s in the arc A . • and f(s,s') be the value of 

the flow s from N„ to N_ , . One of the problems in multi- 

s s 

commodity flow is to find, 


XDBX ^ f (s, s' ) 
s=l 


xf . - ^x 




J 


(5.1.1) 

-f(s,s' 3 if g’ = s, 

0 if 3 4 - s, s' (5.1.2) 


\ w ^ 7- : 

f(s,s' ) if 3 = s’ 


s . t . 
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10 


::^o 



(5.1.4) 


When there is only one kind of flow, an undirected arc can 
always he thought of as two directed arcs of opposite direc- 
tion, aiid arc flows of opposite directions can be concolled. 
But arc flows of different commodities can not conceO. each 
other. This is one of the centro.1 difficulties of the multi- 
commodity flow problem. 

As said- before a single commodity flow problem can 
always be cast into a linear program with objective function 
Z = cx subject to the constra-int Ax = b. The matrix A has 
the unimodular propert 3 ’’, and tho optime.l solution has alwsiys 
integer components. This is no longer true in the multi- 
commodity flow problem. Most multicommodity flow problems 
cannot be solved b 3 r the labelling method or its variations. 

5.2 Mlimm COST Min, TIOOI-mODITT ILO’-J PROBLEMS 

The minimal cost multicommodity flow problem con be 
formulated in both node - arc and arc -- chain form. The 
arc - chain formulation for minimal cost multicommodity 
flow problem from Tomlin £54] is presontod bolow. 

In the formulation it is assumed that there is only 
one source and one sink for each commodity. 
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Ilomencl0.ture 


D . 

1 

k ^k 


arc i in the network 
k 


r — i q 2 B • » • 4 111 « 


• • B c 


ITk 


chains of arcs joining source and 


sink for coiniaodity k. 
Capacity of arc a^« 

Cost per unit flow on s-rc a^. . 


k 

r 


^k 

Q/ • • 
10 


k 


flow of commodity k in chain c'^ 


1 if c, 

0 otherwise 


Vj — lf»BB; ) 


V 


A" 


0. 


m 

Formulation: 


[a!j''-j] arC“Chain incidence matrix, 
i = 1 , . . . v m 

j = Ij.BBy k 


Total number of commodities, 

Requirement of commodity k. 
A very large number b 


Objective Function: 

Objective function will be minimization of total 

cost TC 

ki. lu i_ I- 

min TC = X 2 2. o x'? (5.2.1) 

k=a. j=l 1=1 ■' 

Constraints: 


Capacity constraint 


N- 

2 

k=l j=l 


i .k ,k 

# • • 

10 0 


ar. X- ^ b. 


(i — ly.BBy m) 


( 5 . 2 . 2 ) 



m 
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Flow requirement constraints 

% \ 

Zl X = r, Ck = 1,..., q) (5.2.3) 

j=l ^ ^ 

Tlio main objective of presenting above formulation 
is to illustrate tlie largeness of the linear programs for 
multicommodity flow problems. Norms,lly the special structure 
of these programs is utilised in a.lgoritlims for their 
solution. [ Gremeans, Smith and Tyndall [9]? So.igal [48]] 
Hartman and Las don [24], Kennington and Shalaby [34]. 

How before we proceed to the capacity exp<ansion 
problem for multicommodity flow networks we doscribe the 
multi cominoditj.?' flow feasibility problem (as described in 
Hu [ 30] ) . 

5.3 IfULTICOMIODITY FLOW FEASI3ILITT PRO BlEIIs 

In the feasibility problem of miiLticommodity flows, 
the flow value of each of the commodities is prescribed, 
and the question is whether or not all these flow values can 
bo realized simultaneously in a given network. For example, 
let the network be as shown in Fig. (5.3.1) and the four flow 
requirements bo as shown in Fig. (5.3.2) (i.o. , we need two 
units of flow from to H 2 , three units of flow from I 2 

etc.). ¥e try to ship the four commodities in such a way 
that no arc capacity of any arc in the original network is 


exceeded 
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AsBume that we te-ve a liat of feasible networks that 
satisfy the flow reQ.uirements, These feasible networks all 
have the sa.ine number of arcs connecting the sane pairs of 
nodes as the original given network, but every network has 
different arc capacities, (Some arc capacities may be zero, 
but we still count them a.s arcs of that network, so ths-t all 
feasible networks have the same numb cmr of arcs.) These 
feasible networks are constructed as follows. For each flow 
requirement we find a, chain from the source to the sink of 
that commodity. Then we assign to each arc of the chain the 
capacity equal to the flow requirement. Superimposing all 
the chains for all flow requirements, we have a feasible 
network. . ■ 

There are millions of such feasible networks, each of which 
satisfies the proscribed flow requirements. Any linear convex 
combination of these networks also satisfies the prescribed 
flow requirements. The approach described below determines 
if the original net work contains a network which is a linear 
convex combination of these fea,siblo networks. 

Let [a- • ] be a matrix of m rows, each row correspond- 
ing to an arc of the given network. (In the case of Fig. (5.5.1) 
the matrix would be of five rows.) A column j of the matrix 
represents a network which has as the capacity of the 

first arc, b . 2 ^ as the capacity of the second arc, etc. Then 
all the feasible networks mentioned earlier can be represented 




by columns of the matrix. If the original network represen- 
ted by [b ^5 ...» b^J contains a network [b^,..., b^] with 
b^ bi [i = 1 ,..,, n], and [b^,..., bj^] is a linear convex 


combination of the columns of then the original 

network is feasible. Let the coefficient of the convex 
combination of the column 3 be x.. Then the feasibility 
problem can bo formulated as 


max © 


subject to a^^- x^ + s^ = "*^1 = 1 ,..., mi 3 = 1 ,... 


Xj 0, 0- 


(5.3.1) 


If the optimum value of 9 is 1 , then the original 
network is feasible. There are millions of columns of 


columns in [ a^^j] in this formulationi fortunately, we i 
not write down 0 .II the columns, as the technitiues like 
revised simplex method will prove to bo of much help. 


we need 
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5.4 miHEIIATICAL PROG-EiVM-mia FOEMUIATIOITr 

In this section we present the matheniatical progra.inEiing 
f ormul8,tion for the multicomraoditj’- flow network capacity 
expansion prohlcn in arc-chain form. 


O h.iective I~ nnc tio_n ; 

The objective function is to raiiiimize the total cost 
of introducing va,rious a.rcs at various levels into the 
network. Cost of introducing the k-th level of a^^ can be 
expressed as d^j^ is 0-1 variable. The total cost 

function TC, can be expressed as. 


TC = ^ ‘^‘il ^il '(5.4.1) 

Objective function (5.4.1) is to be optimized considering 
the following constreunts . 


Constraints 


1.’- 


X 

j=l ^ 


= r 


k=l j=l 


k 


10 0 


(li ~ ol) 


^11 yi 1 “ 


il 


(5.4.2) 

(5.4.5) 


The constraints (5.4.2) are flow requirement 
constraints where as constraints (5.4.3) ensure that value 
of y.T is 1 when sum of all commodities flown over a^^ 
exceeds of y ^2 -i-S 0 otherwise. 
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5.5 BMUCH AND BOUITD I^ETHOD; 

/The problem (5.4.1) - (5.4.3) can be solved by any 
of the existing 0-1 programming algorithms. 

The branch and bound method which we have developed 
is exactly same as that of for single commodity flow 
problem, discussed in Section 4.3 except that in Stop 5 
we solve the multicommodity feasibility problem (5.3.1) 
with various arc capacities and commodity roq.uirencnt vector 
as the input, in place of finding maximal flows. 
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C^Ll TIFECPCJ 
CC 220 I»J,MK 

iiiv(ii«ic(n 

22C CCMlKuE 
CC TC 242 

240 ICfPKJ-ICCl'KHl 
1F(IC{RKJ,IE,NCJ€0 TC 2oS 

242 MdCK*! 

241 PPI^E»PK-KE4CK 
IFIPPKEJ5C0tSCCt247 

247 IctPPKE l-lCfHKRE HI 

IF(lCIRKKBl,CT*NiJGO 7C 245 
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CC TC 502 
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540 CCRTI8UI 
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,.cm TiPBiRwi '■■■■• 
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IME6ER C/lC20.20),E(K)!;ilJif;?L , 1,1 
IMECEP CD(W.iSj ®”'’****^*'^*^°*^®^ 
Jisi-q 
EC 1 l*lfK 
EC 3 J««3tN 
CCIiJI-CCCl, • J 
l»l ^ ', ' 

CC 1C 1-1,10 
ICI)-0 

CC 30 4-lfiq 
CKltJl-O 
CCNlIhUl 
K-1 
I«l 

CC ICO 4-ltK 

lE(CClt4)«EC«0}Go t(j |CQ 
IC43«C(I,41 
BUl-l 
CiitltKl-J 

IF{4-EC.MICC 10 3C0 
P-K41 
CCNUNUi 
l-Ml 

irii,ci,W42jcc TO 20 

usu-c 

PC««.i 
f(»l 

I-Cllt-l.F'CM 

IF(I.EC* 0 » 6 C tc 250 

EC 2CQ d«lt^ 

IFCCfl.JJ.ll.OlOC TO 158 
IpdtJUNC.OlCC TC 2CC 
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C4EI Pl ^ iC j 
E(4l-I 
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C <t<1 fit 1*E i tn.KT 1-*# (NT I 
C (NT tlil 1 •€ (NT ,t it U f i NT I 

CtPitiiAT l«C(NW,LI1>|(Nt I 


lOO 

250 


10 1 


158' 

200 


500 


305 


t44(NT } 

iFCf'i'.EC.NSHO TC ?9 

€C 1C 3C5 
2C PPim9tJ^Slfc 
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E^C ^ " 

IIEF1C PIM 

StlEPCuTINE F^n4(KATf |(4TtmINilC} 
iFCPilt.LT.HniCO TO Ic 
HF4C-FJ1T 
CC TC 2C 
IC FI^iC■KIT 
2C PCTUPft 
fUE 
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MAXIHUN* NIImPFR OF ARCS WHICH CAN eR CCNFinFRFr IS 50 
KAXliVlJM KO'. nr i FVELS FOR EACH ARC. IS 4 


nTHFNSTTK STfR A f ^0 ) , j A ( 5 0 ) , 1 8 ( SO ) . IC ( 5 0 ) 
r.TPFKSICN C(?n. 90) .K1 (20,20), rP(20.20> 

Dl KFkS ICN NX (9n, 90 ,1 f!) 

RFAl JljHP( 20 .?n. 10 ) 

INtFgFR R(20.9fi) .FnLi (20,2 o) ,X(20.20 ).cAp(?0,9n) 
INTFRFR tp.rfofi ,sukx 
IMFGFR PRn 
nnST=n.o 
NY7s ?0 

5?n F0R!^AT(3l2*F6'.9.f4. Al ) 

fiOO FORf^ AT (20X,3iA) 

RFAp9oO,DaSh,nPR 
90 0 FORMATcAi, 13) 

Of) afl KKKaS,,NPR 
Re4Cai 0, REoFi,. fi. nS .N f, PRN 
PRTNTaOR.REQFl .N.NS.NT,PRN 
80 9 format (1 H ,3 014) 

51 0 FnRA*AT(aoi4), _ 

R! Gain 00 00 ’.0 

Ksl 

no i 55 I Jjl ,50 
• STFpA{n=o '.0 
15 5 CONTINUE 

no 40 I si, N . 

RFAt^aO, (FOLin* U) *N ) 

FRTnT625. (Foii (j. J) .J=1,N) 

625 FORNaT (4 {Ip ,3nT 4) ) 

690 Fork At ( 20 i 2 ) 

DP 40 jri.N . 

CU, J) *B IG 
CP (I .J)*C( I « j) 

X( I » J) *0 

NX( !« J,K)sinxfRTG) 

JUKP( I, J,K)=RTr 
Kl ( I, J) »1 

CAP( I. J) SB IG 
4(1 CO NT IN UE 

7 ReA!:52 0. I. J.K. JUMP (T .J,K ),NX (I ),SlG 

pR|NT630,I I IM P(|,J,K),NX{I*J»K)»SlG 

630 FoRHAtClH^ ,3T2.F<^'.2. M,Ai) • ' / 

iF (S iG'.eq. dash )R n TO 3 ' ■ ^ 

Go fO 2 

3 COnT INUE 
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rr 11 "! Kai,?h 

If (Rf^FI .GT.FXfT .J.KUGO To g 
?1 K1 ( T . J ) sK 

rr F k^bI ,K 

5 Sji pFsSji;f5iF+ J ump ( j .J .k2) 

6 rn.:.)sSjUMP/Frn4TfRFQFL) 

CF (T .J) = C( I. J1 

RjljM p= n. 0 

Qf: TC 2 q 
1 n rOhT jNiJf- 
F lFfK.F 0 . 4 )Go fn ?1 
2n CHM tkde 
3 0 CnMTNliF 

CAM ShHrT (C .My7 .R 5 

Is1 

TPrNR 
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K>'^*Ki f Jp .NtPi 1 

G0ST=CGST + C( Tp,KiTp.l ) «FL0 At (REoFl ) 

STfPa( DsJUMpf TP.KTF1 ,K3) 
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iBcDsNTPl 
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. PPfK't2oO0,COST, .. . 
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MflSsl 
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suf#x»suHx*xn.NT) 

410 CnNfiNUB , .. 
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3?p m NT In up : : : 

3H P CHlvT iNijE ■ ■■■ ■ 

3q:i PR INT 6t 0, COST.. , , . - 
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C SuRqnuTiJ^P c:HnRf FINDS SHflRTFST CHA iNS BFTWFFN At I PA! 

c ^n^FS IN A nftwfrk 

c 

c 

Fi]nRni;T[^F SHOPTe n.i .nyZ.R) 
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NATN program for, for HEURlRtlC 3 ( SF CT In N . 4 . :* ) 

MAXIMUM NIIMRFR OF NODES WHfCH CAN BE CONSinFRFf JS 10 
Ma !< I Mljy^ MllMRFR QF ARCS WHICH CAN BF CONRIDFRFr TR 50 

Maximum no. of ifvEls for each arc is a 


c 

c 

DtMfNrIoN STEpAf'inl, Ia( 50) 'iBi^O). IcCsO) 
nlyEisSiO^! Cfpf’ .Pnl .Ki (20,2o) ,Cp{20,2Q) ■ 

dtmfnstqn Nx(2n.?n.in) 
real juMp( 20 .?n. ?n ) 

INTfGFR R(:2o#pn^) .FOj 1 (2o.20) ,y ( 20 . 20 ), CAP(2o.?n) 
INfEGER TP.RfOFI .RllMX 
InTpGFR prn 
crsTsO .0 

f^X7s20 , , . , 

■ 520 F0RMAT(3I2,F6.P. f4.Al ) 

600 F0RHAT{20X.3I6) 

REAOrOO.DASH.NPR 
90 0 FORM AT (Al, 13 ) . . 
rO BP KKK*1*NPR 
PFAD51 n,RFOFL..N,NR.NT,PRN 
pRTNT8n5 .RpOFI . N' . NS. N T , PRN 
8(19 FoRyAiIlH ,30 t4) 

51 0 FCrM aT (2nU) 

RiGsinoooo.n 

Ksi 

Do 155 151.50 
STfPaI DsO'.O 

155 CnMlNUF 

on 4(1 T 5 1, N 

RFAn620, (FOIL/ l^. 1.1 ,jsi,N) 

PRtNT 6 25, (FOLl. (I.. J ).>!•» N) 

625 Format < 4 flH »30 t 4)1 

62n FORM aT,{2oI2) 
no 40 vJsl.N 

cc I. j) sBiG , ; 

Cp( I . J) sC ( I . J 1 ' 

xf T . J) sO : - ■ 

K X ( I , J , K 1 5 T F i X ( B f G 1 
JllMpC .1* U» K) aRlO 
K1 ( I . J) =1 
C*P< I* sBtG 

40 CnNtlNliE - 

2 RF 40520. I. J. K. .liIMP <T . J,K ), NX (T .J.K ), SIG 
FR!NT630. I . j. K. jOMPf i ,J,K)»NX(I.J,K).SIG 
630 fORMAT(1H .3I2.F6.2. t4,A1) 

IF f S TG .EO'.DaSHIGO T(1 3 


Gn Tf! ? 

3 CnMTNUp 
SJiJf'psn. 0 
nn 3f| 

nn_2n j=i,n 

nr 10 Ksi . 2 n 

TFfnpQFL.nr.K^fT .j.i<5)G0 To 9 
21 Ki'fT.J^sK 
nn K K?si .K 

5 FjIIN'PsSjL'MP+ J|'|MP( T , J ,k2) 

6 Cf T . j)=SjUMP/Fi nAT(RFOFL) 

CP<T.J)=CCT.J) 

S‘Jii^p=n.n . 

Gn TC 26 

9 IF (K ,Pt3. 4)00 fn pi 

10 Col^TlN^Jp 

26 TF (FpL|, C I. j)'.Fo. 1 ) xn,J) = REQF| 

?0 CtlKflMUe 
3fl cHfgT iK'ue 

CAl L IM (X.N ) 

CA| L SH0RT{C.N.t^Y7.R) 

I»i. 

TP?i^S 

31 MPlsRlfP.Nf) 

K3rKl (TP.NfPi) 

CnSfsCnST + Ccfp.NTpi )«pLOAT(REciFl ) 

IF f K3-1 ) 44.44*46 . 

44 R1FPA( I)=JUMP(Tp.^!Tpl,K3)/FL0AT(X(Tp,NTpl) ) 
rO Tn 48 

46 PTFpA(T)=jUYPfTp.MTp1,K3)/F|_0AT(X(tp,|gTpl)-^Xi’Tp,^iTp1,K3-1)) 
4R I a fl )sTP 
IR(T)sNtPi 
lC(n = K3 
1*1 + 1 
TP#^fP1 

lF(TF.NF.Ivt)Gn fn 31 

lSfUil-1 

PR I^f?OOO.CoSf 

2nOn FnRHAT(iHQ,in<pw--).Fi0‘2) 

CAl.i, DECEND(STFPA . TSTU, lA, iBf fC) 

K0S51 , 

32 Do 5C lal, ISTt'l 
LA*Ia( I) 

LB=TB( I ) , 

LFsTnCT) 

ip ( T .pD. MoS 5 Gn fn 34 

of i A .10) aO'.n ^ 1- V 

CPCI A,l.B) = C{LA..t'R.). 

CAp ( 1, A ,Lp) =K'X fi A ,1 fi. i C ) 

Gn TO 5o ■ ' 
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I f' -t r„ ! ■ ; 

-I- T si ,1. t! , 

■11 Jl|M+!il|MP (I A, lb ,I I) 

'’M A .I.B^-'3iiHj||M>ri.nAT(NX(LA,LB .LC) ) 

f’l" (I \ , I « ) sHd. A ,i'q ) 

HA'Jf I A .1 BJrN'Xd A .i R.I C) 

Kid A.i, msin 
TC (T )=LC 
'5(1 ("OMrAiUP 

1 L n Fi ? r s 1 , fv 
nn A? j=i.N 
C(T,j)=nP(T.J) 

rAii. SwnRKn .k%kY7.r) 

P9n TPaKS 

MP1 5R(TP. NT ) - 

30n IP (FCLI. f -Po.i' )G0 To 30«5 

TFfrOLL (TP,NTP1 1 .FO', (-1) ) GO TO 3i0 
305 XfTP.NTPl IsXcTp.NjPI )+l 

T Fa( TP ,nTP1 1 '.nF*. nAP( TP ,nTP1) )gO to 3 07 
Cf TP.NTPI )=Bin, 

CPCfPfNTPI )sCf Tp.KiTM ) 

307 TPuNTPI 

I F([ TP'.PO, NT )Gn fn 400 
NTPl5R(NfPl .NTl 
on tq 30 0 

31 n )r(^TP1 .TP)sX(Nfpi .TPl-l 

IFfVfNTPt.TPl'.NF.nlGn TO 30? 

C(TP,NTPl)=BTn„_ 
rpaP.NTPl IsCffp.NTPI ) 

GO TO 30 7 
35(i CGNTTNUF 
400 SliPXsO 

nn 4 lO t si ,N 
pliNVsFUNK+X C T .MT 1 

41 n rnNTTMjF , , •' 

fFfSLHX.GF.RFGFI inn TO 360 

or 30 01 : 

36 0 rnoTsn'.o 

POSisHOs+t , 

pn Tf' t4ii ;prn .MnPi 

411 rOPpATCIH ,«PRnpiFW NO./SOLUTThN NO. =«. , T31 
rn 3pn tal.N 
rn 380 j=i,N 

TF f X ( ! . J) . En'.n)Gn^,TD 380 
HRl NT6nO. I . j. x( T , J) 
fO 3 75 K »1 , 1 0 _ 

‘FfK.pQ.XlGO fn 374 

K T s a' ’> 1 

7r (X(T.J)'.LF.NXn,J.KT))GO TO ' 


rr‘?^TsCOST + JlJMP( T . J.K) , 

y ? '1 ''f! \ T i\'in3 , ' ’ . 

"^^Trl\uJF . 

f‘’n;\T 61 n. COST 

(iMa.flTnfAi rr?;Ts«,Fio.2/50(?H. . ) ) 

f'rSsvns + l 

IF (f'FS .GT. iSTlMan TO 7 
" A I I f ^ T U . N ) 

no TO 72 

? CaI I T T^Ff Ni 4 ) ‘ 

fOkT tMIIF • 

s rnp * 

F K n 

c 

G . . 

C StlRROUTTN'F SHORT FINDS SHORTEST CHAIN'S BFTWEFN Ail PAIRS OF 

C NODtS IK A network 

c ■ ■ ^ 

c. 

SIIPROUTINF SHORTin.i ,NYZ»R) 
dihfkSion nip'H.phi 
intfofr R( 2 n, 2 fn 
R I c ? « 1 0 0 0 fi n n 
no tio , 1 * 1 . , 1 . 

no 1l 0 Jal , 1. 

iF(n(i,j)'.eo.pTO)on fo 109 

R(T.J)bJ 
GO TO 110 
1 OR n ( T .vj IsNYZ 

iin oontinuf 

no 140 Jsl.L 
ro 170 T=1,L . 

IFC I .FO. JIGO TO l3n 
on 1 ?n K3i ,!. 

I F( K. FO . J .op', k'. FO . T 1 GO TO 120 
TsD( !» J) + 0( J.K1 
iFcOf I. K)-T)12n. 120.112 
112 Hr t.K)sR(T,J) 

C( I .KisT 
1,20 COKTTNUF 
13 0 COKTTNUF: 

140 fOK Tl NDF 
RFfllRK 

FrvO . . : 

C ■ ■ ■ ; 

C 
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